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The order of the post-Newtonian expansion needed, to extract in a reliable and accurate manner 
the fully general relativistic gravitational wave signal from inspiralling compact binaries, is explored. 
A class of approximate wave forms, called P-approximants, is constructed based on the following 
two inputs: (a) The introduction of two new energy-type and flux- type functions e{v) and f{v), 
respectively, (b) the systematic use of Pade approximation for constructing successive approximants 
of e{v) and f{v). The new P-approximants are not only more effectual (larger overlaps) and more 
faithful (smaller biases) than the standard Taylor approximants, but also converge faster and mono- 
tonically. The presently available O ^-accurate post-Newtonian results can be used to construct 
P-approximate wave forms that provide overlaps with the exact wave form larger than 96.5% im- 
plying that more than 90% of potential events can be detected with the aid of P-approximants as 
opposed to a mere 10-15 % that would be detectable using standard post-Newtonian approximants. 



I. INTRODUCTION AND METHODOLOGY 



Inspiralling compact binaries consisting of neutron stars and/or black holes are among the most promising candidate 
sources for interferometric detectors of gravitational waves such as LIGO and VIRGO. The inspiral wave form enters 
the detector bandwidth during the last few minutes of evolution of the binary. Since the wave form can, in principle, 
be calculated accurately, it should be possible to track the signal phase and hence enhance the signal-to-noise ratio 
by integrating the signal for the time during which the signal lasts in the detector band. This is achieved by filtering 
the detector output with a template which is a copy of the expected signal. Since in general relativity the two-body 
problem has not been solved the exact shape of the binary wave form is not known and experimenters intend to 
use as a template an approximate wave form computed perturbatively with the aid of a post-Newtonian expansion 
p]-pl|. Thus, template wave forms used in detection will be different from the actual signal that may be present in 
the detector output. As a result the overlap of template and signal wave forms would be less than what one would 
expect if they had exactly matched. 

In this paper we explore the order of the post-Newtonian expansion needed to extract in a reliable and accurate 
manner the actual, fully general relativistic signal. Previous attacks on this problem suggested that a 

very high post-Newtonian order (maybe as high as jc' beyond the leading approximation) might be needed for 
a reasonably accurate signal extraction Our conclusions are much more optimistic. We show that, starting 

only from the presently known (t;/c)^-accurate (finite mass) post-Newtonian results [^^10, but using them in a novel 



way, we can construct new template wave forms having overlaps larger than 96.5% with the "exact" wave forms. 
Since a reduction in signal-to-noise ratio by 3% only results in a loss in the number of events by 10%, and since our 
computations indicate that the new templates entail only small biases in the estimate of signal parameters (see Tables 
^ and [X below) , we conclude that presently known post- Newtonian results will be adequate for many years to come. 



Before entering the details of our construction, let us clarify, at the conceptual level, the general methodology of 
this work. Central to our discussion is the following data analysis problem: On the one hand, we have some exact 
gravitational wave form (t\ \h) where \k, k — \, . . . ,n\ are the parameters of the signal (comprising, notably, the 
masses mi and vni of the members of the emitting bina ry ||l6|] ). On the other hand, we have theoretical calculations of 
the motion of p7| , and gravitational radiation from [|6|-|l0|, binary systems of compact bodies (neutron stars or black 
holes). The latter calculations give the post-Newtonian expansions (expansions in powers of u/c) of, essentially [p^ , 
two physically important functions: an energy function E{v) and a gravitational flux function F(v) (see exact defini- 
tions below). Here, the dimensionless argument v is an invariantly defined "velocity" related to the instantaneous 
gravitational wave frequency /^^^ {— twice the orbital frequency) by 
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where m = mi +7712 is the total mass of the binary. Let us denote by Et„ and ^V,, the n* -order Taylor approximants 
of the energy and flux functions, 

n 

ET^=J2Ek{v)v'' =E{v)+0{v"+'), (1.2) 

k=0 



E 

fe=0 



(1.3) 



where 



mi TO2 
(mi + m2)^ 



(1.4) 



is the symmetric mass ratio. For finite 77, the Taylor approximants (L2), ( |l.3| ) are known for n < 5 1 17|, | 6| -|l0| . In 
the test mass limit, 77 ^ 0, E{v) is known exactly and F{v) is known up to the order n ~ 11 [p]-p|ll[. [There 
are logarithmic terms appearing for n > 6 that we shall duly discuss later, but in this Introduction we simplify the 
notation by not introducing them.] 

The problem is to construct a sequence of approximate wave forms h^{t;Xk), starting from the post-Newtonian 
expansions (1.2), (1.3). In formal terms, any such construction defines a map from the set of the Taylor coefficients 
of E and F into the (functional) space of wave forms (see Fig. ^. Up to now, the literature has considered only the 
most standard map, say T, 



{ET„,FT,J^hl{t,X, 



(1.5) 



obtained by inserting the successive Taylor approximants ||20| (1.2), (1.3) into the integral giving the time evolution 
of the gravitational wave phase, see e.g. |0,|l3|. [Details are given below.] In this work, we shall define a new map, 
say "P" , based on a four-stage procedure (Fig. |l]) 



(Et^^FtJ ^ (eT„, /tJ ^ (ep„, /pj ^ (i;[epj, F[ep„ , /pJ) -> h^it, Xk) . 



(1.6) 



The two essential ingredients of our procedure are: (i) the introduction, on theoretical grounds, of two new, 
supposedly more basic and hopefully better behaved, energy-type and flux-type functions, say e{v) and f{v), and (ii) 
the systematic use of Fade approximants (instead of straightforward Taylor expansions) when constructing successive 
approximants of the intermediate functions e{v), f{v). Let us also note that we further differ from previous attacks 
on the problem by using a numerical (discrete) fast Fourier transform to compute the overlaps between the exact and 
approximate wave forms. We find that the previously used analytical stationary phase approximation gives only poor 
estimates of the overlaps (see Table |l|) . 

One of the aims of the present paper is to show that the new sequence of templates h^{t; A) is, in several ways, 
"better" than the standard one h'^{t;X). In this respect, it is convenient to introduce some terminology. We shall 
say that a multi-parameter family of approximate wave forms h^{t;ij.k), k = 1, . . . is an effectual model of some 
exact wave form h-^ {t; Xk)', k = 1, . . . ,n\ (where one allows the number of model parameters to be different from, 
i.e. in practice, strictly smaller than n\) if the overlap, or normalized ambiguity function, between {t\ Xt) and the 
time-translated family h"^{t — t; fj,k), 



A{Xk,fJ.k) 



{h^{t;X),h^{t~T;pi)) 



-.0 ^{h^{t; A), h^it; A)) {h^it; a*), hA{t; a*)) 



(1.7) 



is, after maxirnization on the model parameters /ifc |21[ , larger than some given threshold, e.g. max^^ ^(Afc, /u^) > 0.965 
p2| . [In Eq. (1/7) the scalar product {h, g) denotes the usual Wiener bilinear form involving the noise spectrum Sn{f) 
(see below).] While an effectual model may be a precious tool for the successful detection of a signal, it may do a poor 
job in estimating the values of the signal parameters Xk- We shall then say that a family of approximate wave forms 
h^{t; A^), where the A^ are now supposed to be in correspondence with ( at le ast a subset of) the signal parameters, 

'X 



is a faithful model of h (t; Xk) if the ambiguity function A{Xk^ -^fe )i Eq. (1.7), is maximized for values of the model 



parameters A^ which differ from the exact ones Afc only by acceptably small biases [g3| . A necessary [gj] criterion for 
faithfulness, and one which is very easy to implement in practice, is that the "diagonal" ambiguity ^(Afc,A^ = A^) 
be larger than, say, 0.965. 
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Using this terminology, we shall show in this work that our newly defined map, Eq. ( |l.6| ), defines approximants 
which, for practically all values of n we could t est, are both more effectual (larger overlaps) and more faithful (smaller 
biases) than the standard approximants Eq. (1^). A related property of the approximants defined by Eq. ( |l.6| ), is 
that the convergence of the sequence {h^)nej^ is both faster and much more monotonous than that of the standard 
sequence (ft,J)„gjv'- This will be shown below in the (formal) test mass limit 77 where one knows both the exact 
functions E{v) and (numerically) F{v) [Q, and their Taylor expansions to order v^^ jl^. The convergence will be 
studied b oth "visually" (by plotting successive approximants to E and F) and "metrically" (by using the ambiguity 
function (1.7) to define a distance between normalized wave forms). Most of our convergence tests utilize the rich 
knowledge of the post-Newtonian expansions (p_.2|), (1.3) in the test mass limit i] 0. The very significant qua litative 
and quantitative advantages of the new sequence of approximants, Eq. (1.6), over the standard one, Eq. (1.5), when 
77 ^ 0, make it plausible that the new sequence {h^) will also fare much better in the finite mass case ^ rj < j. 
This question, that we can call the problem of the robustness of our results under the deformations brought by a 
finite value of 77 in the coefficients Ek{ri), Fk{ri) in Eqs. (1.2), ( |1.3| ), is more difficult to investigate, especially because 
one does not know, in this case, the "exact" results for E{v; rj) and F{v; rj). We could, however, check the robustness 
of our construction in two different ways: (i) by studying the "Cauchy criterion" for the convergence of the (short) 
sequence {Hq (rj), /i|'(?7), {rj), (r])) versus that of the corresponding Taylor sequence, and (ii) by introducing a one- 
parameter family of fiducial "exact" functions e;f„ (f), / ^(f) to model the unknown higher-order {n > 6) r/-dependent 
contributions to the post-Newtonian expansions ( |l.2D , (1.3) and by studying for a range of values of the parameter 
Kq the convergence of the short sequence (ft,^(7/), . . . ,h^{ri)) toward the fiducial "exact" wave form h^^{rj). 

Though we believe the work presented below establishes the superiority of the new approximants over the 
standard ones /i^ and shows the practical sufficiency of the presently known u^-accurate post-Newtonian results, we 
still think that it is an important (and challenging) task to improve the (finite mass) post-Newtonian results. Of 
particular importance would be the computation ||25| of the ii^-accurate (equations of motion and) energy function 
in confirming and improving our estimate below of the location of the last stable orbit for 7; ^ 0. Our calculations 
also suggest that knowing E and F to would further improve the effectualness (maximized overlap larger than 
98%) and, more importantly, the faithfulness (diagonal overlap larger than 99.5%) to a level allowing a loss in the 
number of detectable events smaller than 1%, and significantly smaller biases (smaller than 0.5%) in the parameter 
estimations than the present 0{v^) results (about 1 — 5%). 

The rest of this paper is organized as follows: In Sec. ff| we briefly discuss the phasing of restricted post-Newtonian 
gravitational wave forms, wherein corrections are only included to the phase of the wave form and not to the amplitude, 
indicating the way in which energy and flux functions enter the phasing formula. Various forms of energy and flux 
functions are introduced in Sees. Ill and IV, respectively and their performance compared. The ambiguity function, 
which is the overlap integral of two wave forms as a function of their parameters, is discussed in Sec. M and some 
details of its computation by a numerical fast Fouri er tr ansform are given. In Sec.VI we present the results of 
our computations in the test mass case while in Sec. VI] we investigate the robustness of these test mass results 
as completely as possible. Sec. VIII contains our summary and concluding remarks. The paper concludes with two 
appendices. In Appendix |^ we discuss the Fade approximants, their relevant useful properties and list some useful 
formulas used in the computations. In Appendix]^ we discuss carefully the issue of optimizing over the phases and 
provide a clear geometrical picture to implement the procedure. 



II. THE PHASING FORMULA 



To get an accurate expression for the evolving wave form hij{t) emitted by an inspiralling compact binary one 
needs, in principle, to solve two interconnected problems: (i) one must work out (taking into account propagation 
and nonlinear effects) the way the material source generates a gravitational wave, and (ii) one must simultaneously 
work out the evolution of the source (taking into account radiation-reaction effects). The first problem, which in a 
sense deals mainly with the (tensorial) amplitude of the gravitational signal is presently solved to order P-p^. 
Such an approximation on the instantaneous amplitude hij seems quite sufficient in view of the expected sensitivity of 
the LIGO/VIRGO network. On the other hand, the second problem, which determines the evolution of the phase of 
the gravitational signal, is crucial for a successful detection. For simplicity, we shall work here within the "restricted 
wave form" approximation [ p6| , i.e. we shall focus on the main Fourier component of the signal, schematically 
h{t) = a^^(i) cos <^^^(i), where the gravitational wave phase <jP^ is essentially, in the case of a circular binary, 
twice the orbital phase $ : (jP^ {t) = 2^{t). 

We find it conceptually useful to note the analogy between the radio-wave observation of binary pulsars and the 
gravitational-wave observation of a compact binary. High-precision observations of binary pulsars make a crucial use 
of an accurate "timing formula" 
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iPSR 



F[tn■,P^], (2.1) 



linking the rotational phase of the spinning pulsar (stroboscopically observed when (p^^^ = 27rn with n £ A/") to the 
time of arrival tn on Earth of an electromagnetic pulse, and to some parameters pi. Similarly, precise observations of 
an inspiralling compact binary, will need an accurate "phasing formula" , i.e. an accurate mathematical model of the 
continuous evolution of the gravitational wave phase 



,GW 



= 2^ = F[t;p,], (2.2) 



involving a set of parameters {pi} carrying information about the emitting binary system (such as the two masses nii 
and TO2). 

Heuristically relying on a standard energy-balance argument, the time evolution of the orbital phase is determ ined 



by two functions: an energy function E{v), and a flux function F(v). Here the argument v is defined by Eq. (1.1) 
which can be rewritten in terms of the instantaneous orbital angular frequency Q 

V = (ml7)i/3 EE a;i/2 (2.3) 

(as above m = mi + m2 denotes the total mass of the binary). The (dimensionless) energy function E is defined by 

Et,t = mil + E) (2.4) 

where E'tot denotes the total relativistic energy (Bondi mass) of the binary system. The flux function F{v) denotes 
the gravitational luminosity of the system (at the retarded instant where its angular velocity is given by Eq. (p.3[)). 
Note that the three quantities u, E and F are invariantly defined (as global quantities in the instantaneous center 
of mass frame), so that the two functions E{v), F{v) are coordinate-independent constructs. Denoting as above the 
symmetric mass ratio by 77 = 17117122/ {mi + TO2)^, the ene rgy balance equation dEtot/dt = —F gives the following 
parametric representation of the phasing formula Eq. (|2.2| ) (written here for the orbital phase) 

/■"ISO E'{v^ 

t{v) ^t, + m dv -f-^ , (2.5) 



$(«) - + / dvv^ -f( , (2.6) 
Jv F{v) 

where tc and <i>c are integration constants, and where for lisibility we have not introduced a new name (such as w') 
for the dummy integration variable. Note that E'jv) < 0, F{v) > so that both t and $ increase with v. For 
definiteness, we have written the integrals in Eqs. ( ^.5[ ), ( |2.6| ) in terms of a specific reference velocity, chosen here to 
be the velocity corresponding to the last stable circular orbit of the binary. Note that the choice of such a reference 
point is, in fact, entirely arbitrary and a matter of convention as one introduces the two integration constants tc and 
<i>c (which will be optimized later) . The choice Uref = viso is technically and physically natural as it is the value where 
the integrand vanishes (because of E'{v)). The definition (and properties) of our approximants do not depend on this 
choice and the reader is free to use instead his/her favorite reference point. On the other hand, what is not a matter 
of convention is that, in absence of information about the coalescence process, we shall also use viso to define the time 
when the inspiral wave form shuts off. 

The numerical value of viso in the case of a test mass orbiting a black hole (i.e. the limiting case rj — > 0) is i/Ve. 
In the case of binaries of comparable masses (77 ^ 0) v\so is the value of v where E'{v) vanishes. We will discuss below 
ways of estimating v\so{ii)- Knowledge of v\so (considered now has a physical quantity affecting the signal and not 
as a simple reference point) is important in gleaning astrophysical information since the inspiral wave form would 
shut off at that point and the coalescence wave form, whose shape depends on equation of state of stars, etc., would 
begin. One of the questions we address below is whether (as had been suggested |l3|) knowledge of v\so{ii) is crucial 
for getting accurate inspiralling wave templates. 

To warm up, let us recall that in the "Newtonian" approximation (i.e. when using the quadrupole formula for the 
gravitational wave emission) one has 

1 32 

En[v) = --r^v'^ ■ Fn = , (2.7) 

2 5 

so that the above formulas reduce (after redefining the constants of integration, or, equivalently, formally setting 
Viso = 00) to 
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' '^^'^c-^v'^v'^- (2.8) 



The explicit Newtonian phasing formula is obtained by eliminating v and given by 
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$^(t) = - ( - — t ) ; where r = 77^/^ m ("chirp time scale")- (2.9) 



5t 

The corresponding Newtonian gravitational wave amplitude is (for some constant C) 

a%'^iv)^Cv\ (2.10) 

so that the explicit Newtonian templates read 

hNit) = C'it, ~ t)-i cos 1^2$, - 2 (^-^^ • (2.11) 

The crucial issue for working beyond the Newtonian approximation is the availability of sufficiently accurate rep- 
resentations for the two functions E'{v) and F{v). In the astrophysically interesting case of two comparable masses 
orbiting around each other neither of the functions E{v) or F(v) is known exactly and thus one must rely on a post- 
Newtonian expansion for both these quantities. The question is how accurate should our knowledge of the 'energy 
function' E{v) and the 'flux function' F{v) be so that we have only an acceptable reduction in the event rate and a 
tolerable bias in the estimation of parameters. Given some approximants of the energy and flux functions (as func- 
tions of v), say Ea{v), Fa{v), and given some fiducial velocity |^ v^^, we shall define a corresponding approximate 
template 

^h^{t;C,tc,^c,m,7j), (2.12) 
by the following parametric representation in terms of v: 

h'^{v) ^Cv^ cos 2$a(w) , (2.13) 

t{v)=t, + m dv^^, (2.14) 



^Aiv) = a>, + / '° dvv^ . (2.15) 



To compute explicitly /i^(t) we numerically invert Eq. ( ^.14 ) to get v = VA{t) and substitute the result in the other 



equations: h^{t) = CV\{t) cos [2$^(V'yi(t))]. Note that we use the Newtonian approximation for the amplitude as a 
function of v. We could use a more refined approximation, such as an effective (main Fourier mode) scalar amplitude 

u'a^ (v) cx Q.^^ f]^"^ cx v^^ f],!"^ {v) . However, our main purpose here being to study the influence of the choice of 
better approximants to the phase evolution on the quality of the overlaps, it is conceptually cleaner to stick to one 
common approximation for the amplitude (considered as a function of our principal independent variable, v). 

The standard approximants for E{v) and F(v) are simply to use their successive Taylor approximants, Eqs. (|l.2|). 



(1.3). Our strategy for constructing new approximants to E{v) and F{v) is going to be two-pronged. On the one 
hand, using the knowledge of these functions in the test-mass limit and general theoretical information about their 
mathematical structure we shall motivate the use of representations of E{v) and F{v) based on other, supposedly 
more basic energy-type and flux-type functions, say e(v) and f{v). On the other hand, we shall construct Pade-type 
approximants, say ep„, /p,, , for the "basic" functions e{v), f{v), instead of straightforward Taylor a pprox i mants . 



We shall then compare the performance of the various p hasin g formulas defined by inserting in Eqs. ( 2.13 )-( p.l5| ) 
either the standard E'^'^ = Et„, F^'^ = Ft„, Eqs. (O), (O), or the new, two-stage approximants E'X™ = E[epJ, 



pncw _ ^[ep„, /p,J. [In all cases, the approximant of the derivative E'j^{v) is just dEA{v)/dv. 
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III. ENERGY FUNCTION 



Let us motivate the introduction of a new energy function e{v) as a more basic object, hopefully better behaved than 
the total relativistic mass-energy Etot, Eq- (2.4), of the binary system. For this, let us consider the limit TO2/mi 0. 
In this test body limit, i.e a, test particle m2 moving in the background of a Schwarzschild black hole of mass mi, the 
total conserved mass-energy of the binary system reads 



(3.1) 



where A;^ is the time-translation Killing vector, and P2 the 4-momentum of the test mass. [The quantity £2 = ^fc/jP2 
is the well-known conserved relativistic energy of a test particle moving in a stationary background.] At infinity 
— Pi/mi, so that the formal expression of i?tot is -Etot = mi — (pi ■ P2)/mi. This expression is clearly very 
asymmetric in the labels 1 and 2 and has bad analytical properties as a function of mi. Both problems are cured by 
working instead with the standard Mandelstam variable s = E^^^ — — (pi +^2)^ = m\ + m\ — 2(pi • p2). Further, it 
is known that, in quantum two-body problems, the symmetric quantity 



'-2 _ 



El 



TO? 



2toiTO2 



2toiTO2 



(3.2) 



is the best energy function to consider when trying to extend one-body-in-external-field results to two-body results 
p9| . In the limit TO2 <C mi the quantity e reduces simply to e = • P2) / mim2 = £2/^2 + 0{r]). 

In the case of a test mass in circular orbit around a Schwar zsch ild black hole the explicit expression of the quantity 
e in terms of the invariant argument x ^ v'^ = (mfl)^/^, Eq. (2^), is 



1 - 2x 
VI - 3a;" 



(3.3) 



The explicit test-mass result (3.3) suggests that the (unknown) exact two-body function e{x) will have also some 
~ {x — xo)~^^^ singularity in the complex x-plane. This lead us finally to consider, instead of the function e, its square 
or, equivalently, the new energy function 



e{x) = e^ -1 = 



^tot 



2\ 2 

7712 



2mi77l2 



- 1 



(3.4) 



Note that we assume here that the total instantaneous relativistic energy of a binary system (in the center of mass 
frame) can be defined as a time-symmetric functional of positions and velocities (so that E{v) depends on v only 
through X = w^), as the quantity discussed in Sec. VII of Ref. [Q. It remains, however, unclear whether such 

a quantity is well defined at very high post-Newtonian orders and whether it is then related to the gravitational wave 

flux by the standard balance equation. 

Summarizing, our proposal is to use as basic (symmetric) energy function the quantity e{x), Eq. (3.4), instead of 
E{x) = (E'tot — TOi — m2)/{mi + TO2). Given any (approximate or fiducially "exact") function e{x), we shall then 
define the corresponding function E{x) (with x = v^) entering the phasing formulas (2) by solving Eq. ( [3.4| ) in terms 
of -Etot = {mi -|- TO2)(1 -t- -E). Exphcitly, this gives 



E{x) = [1 + 277 + C'ix) - 1) 



1/2 



1. 



(3.5) 



The associated f-derivative entering the phasing formula reads 



E'{v) = 2v 



dE{a 



dx 



vrj 



de{x) 



{l + E{x))^l + e{x) dx 



(3.6) 



Having defined our new, basic energy function e{x), it remains to define the approximants of e{x) that we p ropo se to 
use, whe n one knows only the Taylor expansion of E{x). For guidance, let us note that by inserting Eq. (|3.3[ ) into 
Eq. (3.4) one gets the following exact expression for the test-mass limit of the function e{x) 



e{x; 7] — 0) — ~x 



l-Ax 
l~3x 



-x{l — x — 2>x^ 



•)■ 



(3.7) 
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The generalization of the expansion Eq. (3.7) to non zero values of ry is only known to second post-Newtonian 
(2PN) accuracy. Using Eq. (4.25) of Ref. [§, that is, 
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(9 + 77)x- 



1 / 

27 - 1977 + ^ 



we compute the 2PN expansion of the function e{x) for a finite 77: 



l-(l + |)2:-(3-^77) 



35 



(3.8) 



(3.9) 



The basic idea behind our proposal is that on the grounds of mathematical continuity | |3l[ | between the case 77 ^ and 
the case of finite 77 one can plausibly expect the exact function e{x) to be meromorphically extendable in at least part 
of the complex plane and to admit a simple pole singularity on the real axis oc (a; — a^poic)"^ as nearest singularity in 
the complex cc-plane. We do not know the location of this singularity when 77 ^ 0, but Pade approximants are excellent 
tools for giving accurate representations of functions having such pole singularities. For example, if we knew only the 
2PN-accurate (i.e. 0{v'^)) expansion of the test-mass energy function (6), namely e2PN(a;; 77 = 0) = — x(l — x — 3x^), 
its corresponding 7;"'-accurate diagonal Pade approximant would be uniquely defined (see Appendix 0) as 



ep^{x]r] = 0) = -x 



l-Ax 
1 -3a; 



(3.10) 



which coincides with the exact result, Eq. (3/7). Having reconstructed the exact function e{x), we have also recon- 
structed, using only the informa tion conta ined in the 2PN-accurate expansion, the existence and location of a last 
stable orbit. Indeed, using Eqs. dU) and ( ^ ) we find 



-rjv 



1-67;^ 



(l_3„2)3/2 



(3.11) 



which is the exact test mass expression exhibiting a last stable orbit at 7;iso = 1 / \/6 . In Table | we have compared 

7;^lso predicted by the standard post-Newtonian series and the Pade 



0) and e{x; 77) with < 77 < 



at different post-Newtonian orders the x\, 
approximation to the same. 

It is important to note that our assumption of structural stability between e{x;ri 
is internally consistent in the sense that the coefficients of x and in the square brackets of Eq. (3.9) fractionally 
change, when rj is turned on, only by rather small amounts: 77/3 < ~ -1-8.3% and —3577/36 < —35/144 ~ —24.3%, 
respectively. This contrasts with other attempts to consider 77 as a perturbation parameter, such as Ref. |32|. Indeed, 
in the quantities considered in the latter work several of the 2PN terms have coefficients that vary by very large 
fractional amounts as f] is turned on: some examples being 12 -I- 297/, 2 + 25?; + 27/^, 4 + 417/ -I- 87;^ in Eqs. (2.2) of 
the second reference in [ |3^ . Moreover, the fact that many of the coefficients in their Eqs. (2.2) increase when 7/ is 
turned on (like the ones quoted above) is not a good sign for the reliability of their approach as it means, roughly, 
that the radius of convergence of the particular series they consider tends to decrease as 7/ is turned on. We shall 
attempt below to further test the robustness of our proposal. 

In summary, our proposal is the following: Given some usual Taylor approximant to the normal energy function, 
Et2„ = ~hvx {1 + El x + E2 x"^ + ■ ■ ■ + En a;"), one first computes the corresponding Taylor approximant for the e 
function [B3|, say 



6X2 

in which the only known coefficients are 

ao = 1; «! = 



fc=0 



02 = -3 



3577 

12 



(3.12) 



(3.13) 



Then, one defines the improved approximant corresponding to Eq. ( |3.12 ) by taking the diagonal (P^-, if 77 = 2777.) or 
subdiagonal {Pm+i^ if n = 2m + 1) Pade approximant of —x^^ 6^2,, (x): 



'XP 



^flfc x'' 
.k=0 



(3.14) 
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where e = or 1 depending on whether n = 2m + e is even or odd. For completeness, we recall the definition 
and basic properties of Fade approximant in Appendix Let us only mention here that the P^+e approximants 
are conveniently obtained as a continued fraction. For instance, the Fade approximant of the 2PA^-approximate 
BTiix) — —x{ao + aix + a2X^) is: 



-XCq 



1 + C2X 
.A 



-CqX (1 + C2X) 
1 + {CI+C2)x' 



(3.15) 



By de manding that this agrees with to order v we can relate the c„'s in the above equation to the a„'s in Eq. 
(Ill): 



ai 02 , ai 
Co = ao; ci = ; and C2 = 1 . 



(3.16) 



Explicitly, this gives, 



Co = 1 ; ci = 1 



C2 = -ci 



O 12 V 



4-h + -9v' 



(3.17) 



so that 



(x) = ~x 



1 + i r; - (4 - 



1 „2 



i + i^-(3-i^)' 



(3.18) 



Given a continued fraction approximant ep„ (x) of the truncated Taylor series eT„ of the energy function e{x) the 
corresponding E(x) and E' (x) functions are obtained using: 



Ey>Sx)= l + 2ri[^Jl + evSx)-lj 



1/2 



- 1, 



(3.19) 



2v 



Thus, for instance 



dEp^ (x) 
dx 



VT] 



dep^ (x) 
dx 



(l + £;p„(x)) v/l + ep„(x) 



co(l + 2c2t;^ + (ciC2 + 4)v'^) 

(1 + (ci + C2)«2)2 [1 + Ep,{v^)] Vl + ep,(«2) ■ 



(3.20) 



(3.21) 



where Ep^ is given by Eq. (3.19). The hatted notation introduced in the left-hand side of Eq. (3.21) will again be used 
below and indicat es th at one is dividing some function of v by its Newtonian approximation: e.g. E'{v) = E' (v)/ E'j^{v) 

where, from Eq. (2/7) E'j^{v) = —77 v. 

Having argued that ep^{x), Eq. ( 3.1§| ), and the corresponding Ep^{x) defined by Eq. ( |3.19 ), are better estimates of 
the finite-mass energy functions than their straightforward post-Newtonian approximations, Eqs. (3. J), (3.9), we can 
use our results so far to estimate both the location of the last unstable circular orbit (light ring) and that of the last 
stable circular orbit. The functions ep^{v), Ep^{v) are plotted in Fig. || together with eTi{v) and Et^{v), both sets 
for J] = 1/4, and compared with the exact functions e{v) and E'{v) in the 77 = (i.e. test mass) case. We see that 
the rj = 1/4 P- and T-approximants are smooth deformations of their test-mass limits. Note that the variable x = 
is, in the limit 77 — s- 0, equal to m/r in Schwarzschild coordinates and can be used as a smooth radial coordinate. 
If we wished we could also introduce the function Jtot(a;) giving the cc- variation of the total angular momentum. 
It is indeed related to the total energy Etot{x) by the general identity (for circular orbits) dEtot = ^dJtot where 
the circular frequency is given hy mQ = = x^^"^ . The consideration (even without knowing its precise analytical 
form) of the effective potential for general (non circular) orbits Etot — Etatir-, Jtot) in terms of any smooth radial-type 
variable r measuring the distance between the two bodies allows one to see (by smooth deformation from the 77 = 
case) that the minimum of £'tot(a;) (which necessarily coincides with the minimum of Jtot (a;)) defines the last stable 
circular orbit. Indeed, it is the confluence of the one-parameter sequence of minima of £'tot('', Jtot) considered as a 
function of r for fixed Jtot (stable circular orbits) with the one-parameter sequence of maxima of -Etot(''j Jtot) (unstable 
circular orbits). Note also, from Eq. ( ^.20 ), that the last stable orbit (minimum of i?(a;)) necessarily coincides with 
the minimum of the function e{x). As for the last unstable circular orbit it is clearly defined by the square-root 
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singularity oc {x — Xpoie)~^^^ of E{x), correspondi ng to a simple pole (a; — Xpoic)"^ in e{x). Applying these general 
considerations to our specific 2PA^-Pade proposal ( 3.18| ) one easily finds that we predict the following "locations" (in 
the invariant x variable) for both the light ring (corresponding to r = Sm for a test mass around a Schwarzschild 
black hole), 



''light 1 



and for the last (circular) stable orbit, 



olc(^) = 3 



1 1 



(3.22) 



^Iso 



hi) 



2- 



(3.23) 



We recall that x is invariantly defined in terms of the orbital circular frequency fl — 2tt /°''^ through x = {m fi)^/^, 
so that the gravitational wave frequency (twice the orbital frequency) reads 



/ 



GW 



2/ 



orb 



,3/2 



4397.2 (6x)5 ^ Hz. 

m 



(3.24) 



In the equal mass case (77 = 1/4) Eqs. (|3!22|) , yield Sx^^^^^ (i) = 156/109 ~ 1.4312, 6x1^^ (i) ~ 1.1916, and 



therefore f^^ 



GW _ 



2/j™ — 5719.4(m0/m) Hz. In particular, for a (1.4m0, 1.4m0) neutron star system, we predict 
f\so^ ~ '^fhso ^ 2042.6Hz. Note that our estimate of the (invariant) location of the last stable orbit is significantly 
different from that of Ref. which estimates, for instance — 2f°^'l^ = 1420Hz for a (1.4m0, 1.4m0) system. 
[Actually, we read on the Figures of Ref. || a value (m/°;'')KWW ^ 0.00963 which corresponds to Qx^^"^ ^ 0.925 
and f°J^ ~ 698 Hz (instead of 710Hz marked on their Figures) for the (1.4rn0, 1.4rn0) case.] Qualitatively our 77- 
dependence is different because we find that x^'^^('r]) increases with rj (6a;]^*(?7) > 1 and increasing) while Ref. [3^ ] 
estimates a Qx^^^ {i]) < 1, decreasing with 77. This is an important physical difference as it means, if we are right, that 
binary systems of comparable masses can get closer, orbit faster and emit more gravitational waves before plunging in 
than estimated in Ref. |^2|. As said above, we think that the "hybrid" approximation used in Ref. |^ is not reliable, 
notably because of the stron g 77- dependence (and consequent increase) of the coefficients in their expansion (see also 
the related criticism of Ref. |3J]). We think that our approach (in which the expansion coefficients to e{x) are less 
strongly modified by 77 and where the crucial coefficient 02 decreases with 77 which means a larger radius of convergence) 
is more likely to indicate the correct trend. We have tried in se vera l ways to test the robustness of our conclusions 
under the addition of higher post-Newtonian corrections to Eq. (3.9). We think, however, that such attempts are not 
really conclusive because one does not know in advance what is the "plausible" range of values of 3PN and higher 
77-dependent corrections. [We note in this respect that the range considered in Ref. |32j, [ailmax = |/3i|max = 10, is 
clearly too small as it means, for instance, a fractional change in the coefficient of (m/r)^ when rj changes from 
to 1/4 of 77|q;3|/16 < 16%, while the known fractional change in the coefficient of (m/r)'^ is already 7^29/12 > 60%.] 
In fact, the relative change (ratio afc(l/4 )/Qfc (0) when 77 changes from to 1/4) of the successive coefficients in any 
power series, such as the akijf) in Eq. ( 3.12| ) is expected to increase (or decrease) exponentially with the order k 
due to an 7y-dependent shift of the convergence radius. For instance, in our case if we write the iPN coefficient as 
03(77) = — 9(1 + K3 7/) to model the 'iPN 77-dependence it is not meaningful to consider a priori that K3 can take any 
values in the range ±^2 — ±1 (where we introduced 02(7;) = —3(1 -I- K2 77) with K2 = —35/36). As the negative value 
of K2 has indicated an increase of the radius of convergence with 77 {x^'^y^{r]) = ai{r]) / a2{r]) = i (1 -I- ki 77)/(1 + K2 77) 
with Ki = 1/3) we would rather expect a value of K3 such that 03/02 ^ 02/01, i.e. 1 + K3 77 ~ (1 + ^2 77)^/(1 + ki 77) so 
that K3 ~ —1.9. A value of K3 very different from this estimate (i.e. a value of 03(1/4) very different from —4.8 ) would 
mean that the coefficient 02 (1/4) was accidentally smaller than normal (in which case our estimates (3.22), ( 3.23| ) 
would not be reliable). In conclusion, we think that, given the presently available information, our estimates are more 
internally consistent than previous ones (which include the relevant works quoted in Ref. [^^), but that, if 02(77) is 
only "accidentally" decreased by turning on 77, they might be off the mark. It will be possible to make more precise 
statements on the reliability of Eq. ( ^.23| ) only when the 3PN equations of motion of a binary system are derived 
(or when numerical calculations can reliably locate the last stable orbit). Anyway, we shall see that a knowledge of 
the LSO is not so crucial for extracting the inspiral wave form. [We shall notwithstanding test below the robustness 
of our overall approach under possible uncertainties in the locations of a:poic('7) ^^d a;iso(77).] This is because: (a) 
Interferometer noise rises quadratically beyond a certain frequency; consequently the noise level is pretty high before 
light binaries, such as NS-NS and NS-BH, reach the LSO; only in the case of more massive binaries consisting of black 
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holes and/or supermassive stars with total mass in excess of 25 M©, in the case of initial LIGO, and 60 M©, in the 
case of advanced LIGO, will the frequency at the LSO be in a region where the detector noise is low. In such cases 
it is important to know the location of the LSO accurately because it helps in appropriately truncating the inspiral 
wave form in search templates so that it would not produce anticorrelation with the coalescence wave form which is 
itself not known, as of now, to any accuracy. In the case of lighter mass binaries what is really needed is that the 
approximate energy function should match the exact one at frequencies where the detector noise is the least. This is 
also true for the flux function as we shall sec in the next Section. 



IV. FLUX FUNCTION 

Contrary to the case of the energy function where we could draw on a lot of theoretical information, we have less 
general a priori information on the structure of the flux function F{v). The exact gravitational wave luminosity F 
is not known analytically. It has, however, been computed numerically with good accuracy in the test particle limit 
[J13| and we shall use this in our study. In the test particle limit the flux is also known analytically to a high order in 



(4.1) 



perturbation theory; to order v^^ |11| we have 

F{v;tj = 0) = —v^v^° 
5 

" 11 

J2 AkV^ + {Bqv^ + B^v^ + Bgv"^ + Bwv^° + Bnw") Ini; 
.fc=o 

where the various coefficients Ak and Bk can be read off from pl , 

^0 = 1, 

Ai ^ 0, 

A2 = -3.711309523809524, 

^3 = 12.56637061435917, 

Ai = -4.928461199294533, 
= -38.29283545469344, 

A = 115.7317166756113, 

A^ = -101.5095959597416, 

An = -117.5043907226773, 

^9 = 719.1283422334299, 
Aio = -1216.906991317042, 
All = 958.934970119567, 

Be = -16.3047619047619, 
= 52.74308390022676, 

Bq = -204.8916808741229, 
Sio = 116.6398765941094, 

Bii = 473.6244781742307. (4.2) 

By contrast, in the comparable masses case only the first five Taylor approximants of F{v;r]) are known |]6[-p^. 
Exphcitly, Bki-q) = (A: < 5) and 

^o(^) = 1, 
Ai{v) = 0, 
. / ^ 1247 35 

^^('^)^-W~T2'^' 

Mv) - 47r, 

, , , 44711 9271 65 , 



9072 504 ' 18 
/8191 535 
V 672 ^ '24 



, , , /8191 535 , 

A5(r^) = -[-— + —r^]TT. (4.3) 
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There is, however, a bit of general information about the function F{v) which can be used to motivate the consid- 
eration of a transformed flux function, say f{v), as a better behaved object. Indeed, as pointed out in Ref. Q, the 
function F{v; rj = 0) has a simple pole at the light ring (r = 3m, i.e. x = v'^ = i). The origin of this pole is simple to 
understand physically in a flat spacetime analog. [It is seen from Refs. ||l| and |Q that the curved-spacetime effects 
(metric coefficients. Green function) do not play an essential role and that the origin of the pole can be directly seen 
in the source terms, Eqs. (2.14) of Ref. 0].] Let us consider two (for simplicity identical) mass points, linked by 
a relativistic (Nambu-Goto) string, orbiting around each other on a circle (the string tension T providing the cen- 
tripetal force opposing centrifugal effects) . One can easily find the exact solution of this problem and then estimate 
the linearized gravitational waves emitted by the system ||3^. Let us keep fixed the rest masses nii — m-i — m/2 and 
the radius of the orbit R and increase the tension T so that the particles' velocities v tend to the velocity of light. In 
this limit, one finds that RT ^ p — mw/i/l — ti^/c^ and that the gravitational wave amplitude h cx RT + p ^ p. By 
taking a time derivative and squaring one sees that, as f ^ c, the gravitational flux F ^ fl^ oc p^ tends to infinity 
like (1 — w^/c^)~^. This shows that the finding of Refs. [|l],^ is quite general and that, in particular, it is very plausible 
that a binary system of comparable masses will have a simple pole in F{v) when the bodies tend to the light ring 
orbit. We have seen above that the light ring orbit corresponds to a simple pole Xpoic('7) in the new energy function 
e(x;ri). Let us define the corresponding (invariant) "velocity" Wpoic('7) = \/ x-po\c(ji) ■ This motivates the introduction 
of the following "factored" flux function 



f{v-Tl)= 1- 



F{v;ii). 



(4.4) 



Note that multiplying by 1 — f /wpoic rather than 1 — (u/fpoic)^ bas the advantag e of regularizing the structure of the 
Taylor series of f{v) in introducing a term linear in v (which is absent in Eq. (4.1)). Two further tricks will allow 
us to construct well converging approximants to f{v). First, it is clear (if we think o f v as having the dimension of 
a velocity) that one should normalize the velocity v entering the logarithms in Eq. (4.1) to some relevant velocity 
scale vq. In absence of further information the choice vq = wiso('7) seems justified (the other basic choice uq = Wpoie is 
numerically less desirable as v will never exceed v\so and we wish to minimize the effect of the logarithmic terms). A 
second idea, to reduce the problem to a series amenable to Padeing, is to factorize the logarithms by writing the / 
function in the form 



32 



1 + ln- 



Wlso(?7) 



(4.5) 



The ellipsis in Eq. ( |4.5| ) are meant to represent possible higher powers of In . [Such terms do not show up at order 
v^^ when 77 = and will be also of no concern when considering the ij ^ results at order v^.] The coefficients fk are 
functions of f iso in general [ |36[ | . 

Finally, we define our approximants to the factored flux function f{v) as 



/P„(.;^)^ 1^2.1° 



1 +ln 







' n 




pm 

m+£ 


Ea-' 


\k>6 J 




.k=0 



(4.6) 



where v^^^iv) denotes the LSO velocity (= y/x\so) for the ^"-accurate Fade approximant of e{x), and where P™+g 
denotes as before a diagonal or subdiagonal pade with n = 2m -I- e, e = or 1. The corresponding approximant of the 
flux F{v) is then deflned as 



Fp„(z;;77) = (1 



(4.7) 



where Vp^^^ijj) denotes the pole velocity defined by the ii^-Pade of e(x). For instance, from Eq. (3.22) 



1 / 1 



V5 Vi-i^ 



(4.8) 



Let us now see what this definition gives in practice. In terms of the original expansion coefficients of F{v ), Ak and 
Bk (considered for any rj) and of the fiducial velocity vq = viso, coefficients appearing in the definition (4.6) read 
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^6 


= ^6 , 










= 0, 








h 




A2 Be 








= B9- 


A3 Be 
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= Sio - 


-A2 4 


-Ai 


Be , 


11 


= Bii - 


-A2 £9 


-A3 


4- 



Be . (4.9) 

We find (remarkably?) that in the test particle limit the 0{v^) logarithmic term vanishes identically: 4('7 = 0) = 0. 
The other coefficients are numerically (ry = 0; iipoio = l/\/3; ^'o = ^'iso = 1/^6)' /o — fi — —1.7320508075689, /2 — 
-3.7113095238095, /a = 18.994547272212, ^ = -26.694053570105, /s = -29.756490254383, fe = 196.66395901720, 
/7 = -327.26305863109, /g = 11.063926928123^9 = 1188.0521512280, /m = -2884.9014287843, and /n = 



2823.3603070298. As for the log-factor in Eq. ( [4. 6]) we find that when it is not identically 1 (i.e. when n > 6) 
it is always smaller than about 1.005 for v < viso — 0.40825 and much closer to 1 when v ^ 0.2. Although it is un- 
pleasant to have logarithms mixing with powers, they do not seem to introduce, in the present case (after normalization 
to viso and factorization), a serious obstacle to constructing good approximants to f{v). 

Our primary aim in this work is to compare and contrast the convergence properties of the standard ( "Taylor" ) 
approximants to the phasing formula and its building blocks E{v) and F{v) with the new approximants defined above 
(with their two-stage construction E[ep] and F[fp]). Let us first discuss the case of the flux function which can be 
studied in detail in the limiting case 77 ^ 0. Indeed, in this case one knows both the "exact" (numerical) flux function 
[Q, say Fx{v) and its post-Newtonian expansion up to order v^^ ||ll[]. We can then compare directly the approach 
toward Fx{v), on the one hand, of the succ essive standard Taylor approximants FT„{v]r] — 0) (obtained by keeping 



only the A^. and B^. with fc < n in Eq. (4.1)) and, on the other hand, of the new approximants Fp^ {v; i] = 0) defined 
by Eqs. ( |4.6| ). This comparison of convergence is illustrated in Fig. |3[ We have plotted there, for convenience, the 
"Newton-normalized" flux functions 

FAv)^^^^^V-'v-'"FAv). (4.10) 
Fn[v) 32 

It is clear that the P-approximants converge to the exact values much faster than the Taylor ones. The monotonicity 
of the convergence of the P- approximants is also striking. However, the P- approximants of the flux at certain 
orders (notably v"^ and v^'~') exhibit poles that happen to lie in the region of integration: v\ow < v < v\so- Such 
P-approximants are obviously a bad choice for the construction of templates. Nevertheless, this does not mean that 
one cannot construct P-approximants at that order at all. Recall that in this study we have only considered diagonal 
and subdiagonal Fade approximants of type P™ and Pm+cj respectively. It is perfectly legitimate to employ other 
types of Fade's and in particular the superdiagonal Fade of type P^'^'^- For instance, there is a pole in the region of 
interest in Pj* -approximant of flux while it turns out that the Pg -approximant (which is the one we have used in this 
work instead of Pj* ) does not have a pole in the region of interest. Thus, if one wishes one may trade off a spurious 
zero, in the region of interest, in the denominator of the function with a zero of the numerator, thereby removing 
the troublesome pole (see Appendix ^ for how this may be accomplished via some simple properties of the Fade 
approximants). For completeness we exhibit in Fig. ^ the successive P-approximants to the factored flux function 
f{v;rj = 0). 

The other building block of the phasing formula Eqs. ( 3.10| ) are the approximants to the function E'{v) = dE{v)/dv. 



As we have constructed Ep^ (v) so that it coincides for n > 4 with the exact Ex{v) in the case ry = it would not be 
fair to compare it to the straightforward Et^{v). We need, therefore, to consider the finite mass case rj ^ 0. However, 
in this case, we only know few PN approximations and we do not know the exact result. We can formally bypass this 
problem and have a first test of the robustness of our construction by defining the following fiducial "exact" energy 
function e'^{x): 



"X 



i-(i + !^U- ^'-^^^^^ 



1 — 3(1 ~ Korj) X 



(4.11) 



The 2PiV expansion of e'^{x) coincides by construction with that of the "real" e{x\rf). The parameter enters 
only 3PA'^ and higher order terms. Note that kq parametrizes an infinite number of PN terms in a non-perturbative 
manner because it determines the location of the pole singularity of , namely 
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If we believe our 2PN Pade estimate (3.22) we would expect that a good estimate of the "real" kq (when considering 
77 = i) should be such that 1 - k^V4 = (1 - 35/144)/(l + 1/12), i.e. Kq* = +47/39 ~ +1.2051. To test formally the 
convergence of the sequence of P-approximants away from the region where we know by construction that it would 
converge very fast we shall consider a value of kq substantially different from the Pade-expected one, for instance 
simply kq = which says that the "exact" pole stays when 77 7^ at the test mass value 3x = 1 instead of our result 
above ^.x^l^^ {\) ~ 1.4312. Working again with "Newton- normalized" functions, now 

(4.13) 



we compare in Fig. ^ the convergence of E'j, {v) and E'p {v) toward the fiducial "exact" e'^{v) defined by Eq. ( 4.11 ) 



for rj — J, kq = 0. For completeness we exhibit also in Fig. ^ the successive P-approximants to the "basic" energy 
function e^°. 

The convergence tests performed in this Section have shown at the visual level that the P-approximants behaved 
better than the T-ones. However, the real convergence criterion we are interested in is that defined by overlaps, to 
which we now turn. 



V. AMBIGUITY FUNCTION 



Central to our discussion is the ambiguity function which is a measure of the overlap of two wave forms that may 
differ from each other in not only their parameter values but also in their shape. For instance, one of them could be 
a first post- Newtonian signal corresponding to non-spinning stars parametrized by masses of the two stars and the 
other may be a second post-Newtonian inspiral wave form corresponding to spinning stars parameterized not only by 
the masses of the two stars but also by their spins. Let us therefore consider two wave forms h{t;Xk) and g(t;^k) 
where Xk, k = 1, . . . ,n\, and /itfe, k — 1, . . . , n^, are the parameters of the signals and n\ and are the corresponding 
number of parameters. The scalar product of these two wave forms is defined in Fourier space by 

{Kg){T■\u,^ik)^ J -'^g^Hf;Xkrg*{f;^^k) (5.1) 

where r is the lag of one of the wave forms relative to the other, h{f; Xk) and g{f; ^k) denote the Fourier transforms ]3^ ] 
of h{t; Xk) and g{t]fj,k), respectively, * denotes complex conjugation and 5„(/) is the two-sided noise power spectral 
density. The above scalar product is also the statistics of matched filtering (Wiener filter) which is the strategy used 



in detecting inspiralling binary signals. Sn{f) being a (positive) real, even function of / the scalar product (5.1) 
defines a real bilinear form in h and g. We introduce also the norm \\h\\ = ^ (/i, K). The ambiguity function A is 
defined as the value of the normalized scalar product maximized over the lag parameter r : 

A{Xk,fik) = max ..[^I'^^ii . (5.2) 

^'•^ \M \\g\\ 

where optimization over phases of the signal and the template is symbolically indicated by (f) (see Appendix ^ for 
details). Here Xk can be thought of as the parameters of a signal while Hk those of a template. The signal to 
noise ratio (SNR) for detecting a noise contaminated version of h(t) with a Wiener filter built from g{t — r) reads 
SNR = {h, g)/\\g\\ ■ Its maximum value is SNRmax = {h, h)/\\h\\ = \\h\\ when the time-translated g is perfectly matched 
to the signal: g(t — r) = h(t). Therefore A{Xk,fJ-k) is the reduction in SNR obtained using a template that is not 
necessarily matched on to the signal. 

The dependences of -4(A, /i) on both A and fi are important in designing detection strategies. The dependence on 
the signal parameters A, given some template parameters fi, allow one to define an optimal way of paving the template 
parameter space. The region in the signal parameter space for which a given template obtains SNRs larger than a 
certain value jssj (sometimes called the minimal match ) is the span of that template ji^ and the templates should 
be so chosen that together they span the entire signal parameter space of interest with the least overlap of one other's 
spans. In our case, we are mainly interested in keeping the signal parameters A fixed, and varying the template ones 
/i. In searching for a coalescing binary signal in the output of a detector one maximizes over a given bank of templates 
(i.e. over a dense lattice of /i values). Thus, the quantity of interest is the maximum of the ambiguity function over 
the entire parameter space of templates. This maximum, in the case of identical signals, occurs when the parameters 
of the template and the signal are equal and is equal to 1. However, in reality the template wave forms are not 



13 



identical to the fully general relativistic signal and hence the maximum overlap will in general be less than 1 (Schwarz 
inequality) and would occur not when the parameters are matched but when they are mismatched: 



max^ < 1. 



(5.3) 



If the template wave forms are not 'close' to signal wave forms then it is reasonable to expect that the maximum 
occurs when {Xk — fJ-kl is fractionally rather large. In this case there is not only a substantial reduction in the maximum 
SNR that can be achieved by using such a bank of templates but there would also be a large systematic bias in the 
measurement of parameters. Using the terminology of the Introduction such template wave forms would be neither 
effectual nor faithful. For detection purposes we wish to construct effectual templates, i.e. templates having large 
overlap after maximization over /i. For parameter estimation we further need to construct faithful templates which 
have large overlaps when ~ A. A practical (non rigorous) criterion for faithfulness is that the "diagonal" ambiguity 
function ^(A, A) be close to 1. 

Reduction in the overlap of template wave forms and true signals has an effect on the number of detectable events 
or, equivalently, loss in the detection probability of a signal of a given strength. For a given signal-to-noise ratio, the 
distance up to which a detector can see depends primarily on the amplitude of the wave. Unavailability of a copy 
of the true signal means that the effective strength of the signal reduces from to Ahf) and hence the span of a 
detector reduces by the factor A. The number of events a detector can detect being proportional to the cube of the 
distance, a reduction in the overlap by a factor A means a drop in the number of detectable events, as compared to 
the case where a knowledge of the true wave form was available, by a factor A'^ . For instance, a 10% (20%) loss in 
the overlap would mean a 27% (50%) loss in the number of events |^^. The aim of FN calculations is to make this 
overlap as close to 1 as possible. If we demand that we should be able to detect with FN templates about 90% (99%) 
of the signals that we would detect had we known the general relativistic signal, then we should have the overlap to 
be no less than about 0.965 (0.997). 

As a model for noise we use the expected noise power spectral density in the initial LIGO interferometer pi[: 



Sn{f) 



So 



a + 2 



f>fs 



(5.4) 



where Sq, a and fo are constants that characterize the detector sensitivity, effective bandwidth and the frequency at 
which the detector noise is the lowest, respectively. In the case of initial LIGO a = 2, /s = 40 Hz and fo — 200 Hz. 
Due to the fact that the noise is essentially infinite below the seismic cutoff f^ and sin ce w e terminate the template 
wave forms when the velocity reaches that of the last stable orbit the overlap integral (J^^ reduces to 



A{Xk,f^k)^\\h\\-'\\g\\- 



^ Hf;Xk)9*U;f^k] 



Snif) 



(5.5) 



where f\so is the gravitational wave frequency corresponding to the last stable orbit. In order to compute the maximum 
overlap we proceed in the following manner. The evolution of phase as a function of time is obtained by inverting 
numerically v in terms of t from Eq. (2.14) and inserting the result in Eq. (2.15) and then (2.13). Though, the 



iterative procedure in inverting v in terms of t is rather computationally intensive, we need to employ it since the 
inaccuracies introduced by the stationary phase approximation in computing the Fourier transform of the wave form 
increase with the order of approximation especially in the case of NS-BH and BH-BH binaries. In Table ||, we give a 
measure of the inaccuracies introduced by the stationary phase approximation at various post-Newtonian orders by 
computing the integral in Eq. (5.5) with h{f) being the fast Fourier transform and g{f) being the stationary phase 
approximation (The three cases Aq, Bo and C'o are defined below). 



VI. RESULTS AND DISCUSSION 



Having in hand the ambiguity function to measure the closeness of two wave forms 1 42 we can use it to pursue at 
a quantitative level the analysis of the convergence of the sequence of approximants defined above. 

Let us first consider the wave forms defined in the formal test mass limit where one keeps the rj factors in front of 
E{v) and F{v) but neglects the 77-dependence in the Taylor coefficients of E'{v) and F{v). Explicitly we mean the 
wave forms defined by eliminating (numerically) v between 

h^{v) ^ Cv'^ coa2<i>Aiv) , (6.1) 
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t(v)^U--7j-^m dvV' / , (6.2) 



i>.w=*e-4.-r^..-f^^^, (6.3) 

•j^ -Fa(«; '7 = 0) 

in which v\so ~ wiso(?7 = 0) = 1/V6- Note that the main purpose of the overlap computations made for this formal 
test mass limit is to compare quantitatively the convergence of the P-approximants to that of the T-ones. One should 
keep in mind that when studying below in the formal test mass limit (1.4to0, 1.4toq) or (lOmQ, lOm©) systems (for 
which 1] takes its largest value) the absolute values of the overlaps are not reliable, though one assumes that the 
lessons learned from the P/T comparison are. The absolute values of the overlaps for the (1.4to0,1Otoq) case are 
probably more reliable, but this is not clear as ~ 0.1077 is then only a factor 2.32 smaller that rymax = 0.25. This 
being said we wish to compare semi-maximized overlaps that we can denote for simplicity as 

(r0(mi,m2),X0(mi,m2)) = 

maxf^ <(.x (hl\tc, mi,m2),h^\0, $f , toi, mz)) , (6.4) 



(P°(mi,m2),X"(TOi,m2)) = 

maxj^^^A (tc,^c^mi,m2),h^ (0, $f , mi, TO2)) . 



(6.5) 



Here the superscript on T, P or X denotes the abo ve d efi ned formal 77 = limit of Taylor, Pade-type or eXact 
wave forms, respectively (i.e. A — T,Poi-X in Eqs. (|6.l| ) — (6^)). Here one considers o nly the sam e values for the 
two dynamical parameters of those signals (i.e. the explicit m and 77 appearing in Eqs. (6.1) — (6.3), here expressed 
in terms of mi and m2 in order to psychologically minimize the formal inconsistency of setting r/ = in part of 
the formula and keeping it elsewhere) and maximize over the kinematical ones tf, ^f, , . To maximize over 
the reference times, it is sufficient (as indicated above) to fix i;^ = and maximize over t"^ = tc (t, the time lag). 
Maximizing over the reference phases is more subtle as the overlap depends separately on $^ and and not only 
on their difference. There is, however, a computationally non-intensive way to do it which is based on a conceptually 
simple geometrical formulation of the problem (see Appendix ^). 

Note that in Eqs. (6.4) the approximate template parameters are not optimized, but are taken to be equal to 
that of the exact signal. In other words we compare the faithfulness of the various approximants together with 
their convergence properties. The results are given in Table HI, for n = 4 — 11 43| as well as for the Newtonian 
approximants for the purpose of comparison. The overlaps quoted are the minimax overlaps, Eq. (B12), together 
with the corresponding best overlaps Eq. (Bll) in parenthesis below the minimax overlap. (The P-approximant 



P4 corresponding to n = 7 has a singularity in the region of interest and hence we have used the approximant P3 . 
The Pg -approximant too has a pole and we have not computed the overlaps in this case though if one desires one 
can compute other P-approximants, such as P4 or Pg , at this order.) We consider three prototype cases, say case 
Ao[(1.4mo, 1.4mo)], case Po[(l-4mo, lOm©)], and case Co[(lOm0, lOm©)]. We added an index zero to recall the fact 
that 77 = has been used in E' and F. [One should keep in mind the warning above that the numerical results for 
case Bq are physically more reliable, while Aq and Co are just mathematical ways of testing the convergence.] 

We performed another convergence test (still in the formal rj —^ limit) of a different nature. It is known in 
mathematics that one does not need to know in advance the limit of a sequence to test its convergence. One can 
instead use Cauchy's criterion which says (roughly) that the sequence converges if, given some distance function 
d{h,g), d{hn,hm) ^ as both n and m get large. In our case we have a distance function |4^ defined by the 
ambiguity function and we can compare the Cauchy convergence of the T and P approximants. Some results are 



given in Table IV where one exhibits the semi-maximized (in the sense of Eqs. (|6.4[)) bes t overlaps {T^,T^,^ versus 



(P°, P°+i), for n = 4, . . . , 11, and the three prototype cases Ao, Po, Cq. (As in Table HI where appropriate we have 



used the Pg -approximant instead of Pf. Since the Pg approximant has a pole in the region of interest the entries 
corresponding to n = 10 are blank and the entries corresponding to n = 9 are the overlaps (Pg , Pfi)-) 

The last two Tables show very clearly that the P-approximants converge much better than the T-ones and that 
they provide a much more faithful representation of the signal. To measure the effectualness of our approximants (in 
the technical sense defined above) and study the biases they can introduce, we also performed numerical calculations 
in which we maximized over all parameters, say 



((T°,X°))(mi,m2) 



(6.6) 
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{{P^,X°))imum2) = max„.^„A(pO(mf ,m^),X"(mi,m2)) , (6.7) 

while keeping track of the parameter values mf, which, given the signal parameters mi, TO2, maximize the overlaps. 
The results are presented in Table |^ for the three prototype cases ^o, Bq, Cq and for the most important values (for 
the near future) of the order of approximation: n = 4, 5 and 6. In this case the overlaps are the minimax overlaps. 

Our test mass results sum up the general behavior of the different approximants pretty well. First let us note that 
even at 0(v^^) T -approximants do not achieve the requisite overlap of 0.965 except in the case of light binaries. This 
is consistent with the concern often expressed in the literature about the need for higher order post-Newtonian wave 
forms. In our view the most worrying aspect of the T-approximant is not that it does not obtain a high overlap but 
that the behavior of the approximant is oscillatory in nature. For instance, the 0{v^) T-wave form achieves an overlap, 
with the exact wave form, of about 0.96 which reduces at 0{v^) to as low as 0.71 for system Bq and 0.85 for system 
Co (though for system Aq it maintains a level of 0.965), increases at 0{v^) to about 0.93 for these systems and again 
drops back at 0{v^^) to 0.85 and 0.90 for systems Bq and Co, respectively. One clearly notices that P-approximants 
do not show such an erratic behavior. Recall that, in the test mass case, we are comparing a known exact wave form 
with an approximate signal model and hence the above conclusions are free from any prejudice. Though the second 
post-Newtonian P-wave form is not a faithful signal model, at 5/2 post-Newtonian order the P-approximant is a 
faithful signal model. 



Moreover, P-approximants show an excellent Cauchy convergence as evident from Table IV, Notice that the T- 
approximants have a poor Cauchy convergence for systems Bq and Cq. This makes them ill-suited as faithful templates. 
T-approximants are not always effectual signal models either. Sometimes they do obtain overlaps larger than 96.5% 
but at the cost of producing a very large bias in the estimation of total mass. This is to be contrasted with the 
P-approximants which are effectual at 0{v*) at the level of 99.7% or better at the cost of very little bias {5m /m 
always less than 3.5% and less than 1% in most cases). We have also computed the biases in the estimation of the 
parameter 77 and there too we see a similar trend. 



VII. ROBUSTNESS 



Up to this point in the paper we have mainly relied on the test mass limit to assess the quality of our approximants. 
In this section we shall try to go beyond this formal limit to check the robustness of our proposal under the turning 
on of rj. 

We can first use all the existing information about the comparable masses case and see whether turning on rj modifies 
in any way the trend we saw above. As a first test (a "visual" one) we plot in Fig. ^ the Newton-normalized flux 
functions, Pr„ {v\ rj), Fp^ {v; rj) as a function of v, for the maximal value V — j for the cases where we know them, 
i.e. n = 2, 3, 4 and 5. Using the same information we can also check the 77-robustness of ou r Cauchy-convergence 
criterion. This is done in Table [v| where we present the semi-maximized 6est overlaps Eq. (Bll) {P^ {m,ri), {m,ri)) , 
{P^{m,ri),P^^{m,ri)) and compare them to their T-counterparts for the (real) cases A, B and C. 

We also made many attempts at testing the robustness of our conclusions when taking into account the existence of 
(unknown) higher-order T^-dependent corrections. There is no really conclusive way of achieving such a task but here 
is our best attempt: Our starting point is to model an infinite nu mber o f (unknown) higher-order PN corrections by 
just one (non perturbative) parameter: kq. As introduced in Eq. ( [4.1l| ) above, kq parametrizes our ignorance about 
the true location of the light ring (pole in e{x) and F{v)). Our 2PA^ Fade estimates gave us an ry-corrected value 
^poie) but we wish to consider here the possibility that maybe the true value is quite different from our estimate. More 



precisely Eq. ( 4.11 ) parametrizes the pole at 3a;poic = (1 — x-Qr]) ^, while ^ix^^^^ ~ 1.4312 for = i corresponding 



to Kq ~ -1-1.2051. To explore a very large range of possibilities we shall consider that the true value of kq (for 



rj = 0.25) might range between k,q — ~1 (meaning 3a;poio = 0.8) and kq — +2 (meaning 3a;poic = 2.0). In Table |VII 
we compare the location of the the last stable orbit xiso = w^lso predicted by the T and P-approximants to the energy 
function relative to the exact location x^^. We see that P-approximants capture the location much better than the 
T-approximants. 



Having chosen the range of kq we shall consider, and adopting the definition (4.11) for the corresponding fiducial 
"exact" e-function, it remains to define a corresponding fiducial "exact" /-function, having the property that the 
corresponding P-function coincides, up to 0{v^) terms, with the known expansion of F. To this effect the simplest 
proposal is to define first the Tn (Taylor to v^^) expansion of /^"(w) by 



rn[/^"]= 1- V3(l-«o?7)^ 
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Ak (v) + E (^fc (0) v'^ + Bk (0) v'^ In v) 



.k=0 



k=6 



(7.1) 



where Ak{r]), fc < 5, are given by Eq. 
fiducial "exact" /-function by: 



, and the others (rj — 0) by Eq. Then we define the corresponding 



f^°{v;v) = fp,dv;v) = the R.H.S. of Eq. (4.6). 



(7.2) 



Having defined some fiducial "exact" e and / functions we have correspondingly defined some "exact" wave form h'^ 
and, using the definitions above, both T-type and P-type approximants of this wave form. We are interested in knowing 
whether the P-approximants behave better than the T-ones even in presence of higher-order effects significantly 
differ ent from the behavior expected from the 2PN Pade results. The results of this exercise are presented in Table 
VIII where one has computed the semi-optimized mmimaa; overlaps (P,^(m, 77), X^^^ (m, 77)) and {T^{m,r]), XJl^^m^rj)) 
for the cases A, B, C, for kq 1, 1.2051 and 2 and for n — 4, 5, 6 and 7. In order to test the effectualness 



of the approximants, in Table [X we quote the fully-optimized but minimax overlaps {{P^{'m,r])^ XJl^(m,r]))) and 
{{T^J{m, ?7), (m, 7 7))) again for the cases A, B, C, for kq = —1, 1.2051 and 2 and for n = 4, 5, 6 and 7. 
From Table VIII we clearly see that T-approximants fail to be faithful signal models even at the third post- 



Newtonian order. The second post-Newtonain wave form of this family would clearly fail to capture even 20% of 
all potential NS-NS events that would be detectable with the aid of a family of templates constructed out of P- 
approximants. Even when parameter values are extreme {kq — —1, and very high masses) the presently available 5/2 
post-Newtonian energy and fiux functions are sufficient to construct a faithful P-approximant. 

We observe that except when the parameter values are extreme (very low value of kq and high masses) 0{v^) 
P-approximants are indeed good effectual signal models. In fact in all cases, but one, they obtain an overlap in excess 
of 99%. Bias in the estimation of the total mass is at worst 7.6% and in many cases it is below 2%. On the contrary 
standard second post-Newtonian approximants are not effectual in many cases; when they are effectual they often 
produce a relatively large bias. For instance, for system B, when kq — 47/39, second post-Newtonain T-approximant 
acquires an overlap of 0.98 compared to 1.00 acquired by the P-approximant of the same order. However, the bias is 
97% in the former case as compared to a tiny 1.1% in the latter case. Similarly, for — 2, the 2.5 post-Newtonian 
T-approximant achieves an overlap of 0.988 at a bias of 75% while the P-wave form achieves 0.996 overlap with 
practically no bias at all. The biases in the estimation of the ?7-parameter (not shown) are also pretty small when 
P-approximants are used as compared to T-approximants. 

A word of caution is in order for those who desire to use standard post-Newtonian templates: A careful examination 
of the above Tables reveals that the 2.5 post-Newtonian T-approximant systematically obtains poorer overlaps and 
larger biases. This is of course related to the fact that the 5/2 post-Newtonian fiux is very badly behaved (cf. Figj^). 
Hence one must never employ 2.5 post-Newtonian T-approximant for searches. However, P-approximants do not 
suffer from this predicament. Indeed at 0{v^) P-wave form is an excellent effectual signal models. For all systems 
and parameters this model obtains an overlap of better than 99.5% at a bias less than 1.5%. 



VIII. CONCLUSIONS 



In this work we have studied the convergence properties of various post-Newtonian templates to detect gravitational 
waves emitted by inspiralling compact binaries consisting of neutron stars and/or black holes. We have shown that 
the standard post-Newtonian filters, referred to as the T-approximants that are based on Taylor series, considered in 
the literature define a badly convergent sequence of approximants. Even at order v^^ the T-approximants only provide 
overlaps ~ 0.86 with the exact signal in the case of binaries consisting of 1.4-10 Mq systems. Worse, the convergence 
of the sequence of T-approximants is oscillatory rather than monotonous. Our results on T-approximants confirm 
previous, less convincing arguments in the literature, which were either based on rough quantitative estimates, or on 
numerical calculations based on the stationary phase approximation for Fourier transforms — an approximation that 
we have shown not to be sufficiently accurate for this purpose (see Table ||. 

We have defined a new sequence of approximants, referred to as the P-approximants, based on two ingredients: 
(i) the introduction, on theoretical ground, of two new energy-type and fiux-type functions e(u) and f{v), instead of 
the conventionally used E{v) and F{v) and (ii) the systematic use of Fade approximation for constructing successive 
approximants of e(u) and f{v). The new sequence of P-approximants has been shown to exhibit a systematically better 
convergence behavior than the T-approximants. The overlaps they achieve at a fixed post-Newtonian order are usually 
much higher, and the convergence is essentially monotonous instead of oscillatory (as pictorially described in Fig 
and mathematically measured by the overlaps quoted in Tables HI, VIH, and IX). From our extensive study of the 
formal "test-mass limit" r] = mim2/(mi +7712)^ ^ 0, i.e. keeping overall 77-factors but neglecting 77 in the coefficients 
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of the post-Newtonian expansions, it appears that the presently known 0(w/c)^-accurate post- Newtonian results allow 
one to construct approximants having overlaps larger than 96.5% (overlaps corresponding to kq = 47/39, 2 in Table 
IX and all, but one, overlaps in Table [X) with the exact signals. Such overlaps are enough to guarantee that no 
more than 10% of signals may remain undetected. By contrast (w/c) ^-accurate T-approximants only give overlaps of 
50%, and sometimes even as low as 30%, corresponding to a loss of 87.5% and 97% events, respectively. Our results 
are summarized in Fig. H where we have plotted the fraction of events which the templates constructed out of T 
and P-approximants would detect relative to the total number of events that would have been detectable if we have 
had access to the true signal. We clearly notice the superiority of the P-approximants. Moreover, our computations 
indicate that the new templates entail only acceptably small biases in the estimation of signal parameters (see Tables 
and [X). In the terminology introduced in the text, P-approximants are both more effectual (higher fully maximized 
overlaps), and more faithful (smaller biases) than the usual T-approximants. The above conclusions are primarily 
based on the study of the formal test-mass limit and assumes that turning on rj brings only a smooth deformation 
of what happens at 77 ^ 0. We have also studied the effect of turning on 77 (77 ^ 0) in the coefficients of the post- 
Newtonian expansions. From all our checks it seems that the ry-dependence is indeed smooth and should not alter the 
fact the P-approximants have a better convergence than the T-ones. Our c onstruction predicts that the last stable 
circular orbit is closer (i.e. greater orbital period) when t] ^ (see Eq. (3.23). This is good news because it improves 
the efficiency of P-approximants to be used as filters for detectors having a fixed frequency band. However, we have 
no independent confirmation of this (favorable) dependence on rj. We have tested the robustness of our conclusions 
against possible very drastic changes brought by (still unknown) 77-dependent terms in the higher post-Newtonian 
coefficients. In the case where these extreme changes go in the opposite direction of what is suggested by presently 
known results (i.e. in the case Kq = — 1), we find that the overlaps are worsened compared to our best estimate range 
{kq = 47/39). This shows that it is important to extend the presently available 0{v^) post-Newtonian results to 
the third post-Newtonian level (notably for the equations of motion) . This will allow one to check whether the 
ry-dependence of the 2.5 post-Newtonian results that we use is typical of the higher terms (as our method assumes) or 
exhibit some abnormal behavior for some unforeseeable reason. When third post-Newtonian results are available it is 
still advisable to use the P-approximants: they have consistently higher overlaps and lower biases (cf. see Table IX). 

In this study we have only considered the noise power spectral density corresponding to initial LIGO interferometers. 
Naturally, one must study other cases as well. Based on the current study we can be confident that in all cases the 
P-approximant wave forms will fare much better compared to the standard post-Newtonian ones. However, their 
performance in absolute terms needs to be re-assessed since other interferometers, such as GEO600, VIRGO, and 
enhanced LIGO, happen to have effective bandwidths and the frequency of maximum sensitivity somewhat different 
from initial LIGO. In addition, one must also address the performance of P-approximate wave forms with regard to 
parameter estimation. 



ACKNOWLEDGMENTS 



It is a pleasure to thank Eric Poisson for providing the numerical test mass flux. BRI thanks the Institut des Hautes 
Etudes Scientifiques , University of Wales Cardiff and the Albert Einstein Institute, Potsdam, while BSS thanks the 
Raman Research Institute and Institute des Hautes Etudes Scientifique for hospitality during different phases of this 
work. This work was supported in part by NSF grant PHY-9424337. BSS thanks Kip Thorne and his group for useful 
conversations. 



APPENDIX A: PADE APPROXIMANTS 



A Pade approximant to the truncated Taylor series expansion of a function is a rational polynomial with the same 
number of coefficients as the latter. The coefficients of the Pade approximant are uniquely determined by reexpanding 
the Pade approximant to the same order as the truncated Taylor series and demanding that the two agree. In our study 
we use a continued fraction form of the (near diagonal) Pade approximant instead of the usual rational polynomial. 

Let Sn{v) = ao + ai V + ■ ■ ■ + Un v"" be a truncated Taylor series. A Pade approximant of the function whose Taylor 
approximant to order w" is Sn is defined by two integers m, k such that m + k = n. If r„[- • •] denotes the operation 
of expanding a function in Taylor series and truncating it to accuracy w" (included), the P™ Pade approximant of Sn 
is defined by 

^r(«) = ; TniP^iv)] ^ Sn{v) , (Al) 
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where Nm and Dk are polynomials in v of order m and k respectively. If one assumes that Dk {v) is normahzed so that 



Dk{0) = 1; i.e. Dk{v) = I + qiv + ■ ■ one shows that Fade approximants are uniquely defined by (|Al|). Note that, 
triviaUy, P(5'[S'„] = Sn which indicates that Fade approximants are really useful when k 7^ 0. Actually, it seems that 
in many cases the most useful Fade approximants are the ones near the "diagonal" , m — k, i.e. P™ ii n — 2m is even, 
and Pm^^ or Pm+i n — 2m + 1 is odd. In this work we shall use, except when specified otherwise, the diagonal 
(P™) and the "subdiagonal" {P^+i) approximants. For instance, the P4-approximant of the flux function has a pole 
and therefore we use instead the P3 -approximant. The diagonal (P^) or subdiagonal (P™+i) Fade approximants can 
be conveniently written in a continued fraction form (see e.g. |4q]). For example, given 



one looks for 



and given 



S2{v) ^ ao + aiv + a2V , (A2) 



Pi i^) ^ 1 , cvv = Co 7—7 — - — ^ , A3) 
l + Tife 1+ C1+C2V 



Ssiv) = flo + ai w + 02 + 03 , (A4) 

one looks for 

1 + ° 1 + (Ci + C2 + C3)V + CxC^V^ ^ ' 

The main advantage of using the continued fraction representations is that the lower order coefficients Ck remain 
unchanged as we increase th e or der of the polynomial being approximate d. T his is not true for the coefficients of the 



polynomials -/V„i, Dfc in Eq. (Al). [This is easily seen by comparing Eqs. (A3) and (A5).] The 's are algorithmically 



obtainable in terms of the coefficients in Sn with t < k. For instance, 

Co = ao , 
ai 



Cl 



C2 



C3 



C4 



C5 



Co 

Cocf - 02 
CqCi 

C0Cl(c2 + Ci)^ + Q3 

C0C1C2 

CoCl(c2 + Ci)3 + CoCiC2C3(c3 + 2C2 + 2Ci) - Oj 

CQC1C2C3 

CoCl((c2 + Ci)2 + C2C3)2 + CoCiC2C3(c4 + C3 + C2 + Ci)^ + 05 



C0C1C2C3C4 

C6 = {coCl(ci + 02)^ + 

C0C1C2C3C4C5 

C0ClC2C3[(ci +C2+ Cs)^ + 3(ci + 02)^ + 3c2C3(ci + C2) + C3(c2 ~ Cl)] + 
C0ClC2C3C4[(ci + C2 + C3 + €4)^ + 2(ci + C2 + C3)^ + C3(c4 - 2ci)] + 

C0C1C2C3C4C5[C5 + 2(c4 + C3 + C2 + Cl)] - Oe} (A6) 

The expressions for the later coefficients are long and not listed here. Explicitly in terms of the a's some of the above 
coefficients are , 



Co 


= ao , 






_ Ol 




Cl 








oq ' 






_ 02 




C2 








ai 


ao 



Qo(aia3 - 02) , . „x 

C3 = —r^ r ■ (A7) 

ai(ai — a2ao) 
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A few other properties of the Pade approximants are useful to notice, such as 

pnTnm]^iPL[T4f-']])-\ 



(A8) 



P™+n5„]=ao + ait'P™+J5„-i] 



(A9) 



where (5 = or 1 and where Sn~i is defined by "factoring after w": Sn = ao + oi v Sn~i- Eq. (A9) shows how to obtain 
the superdiagonal Pade of Sn from the continued fraction approximants (= P,' 

Pn'^'^ etc. . . [An alternative way would be, from Eq. (A8), to work with the inverse of the series Sn, 



of Sn-i- This can be iterated to 



APPENDIX B: OPTIMIZING OVER THE PHASES 



The "eXact" (label X) and approximate (label A) template wave forms have the form 

h"" {t; t^,<l>^,C'') = C'' a{t - tf) cos[0f + 0^ (< - tf)] , 

h^{t- t^cbt C^) = a{t ~ tf ) cos[0^ + (l)^{t'~ t^)] , (Bl) 

where we denoted (j)^ = 2$^, (j)^ {t) = 2<i>"^(i), etc. The normalized overlap between and depends on the time 
difference tf — t^ , and (separately) on the two phases </)^, Here, we show how, for any given time lag t — — t^ 
(i.e. after having fixed, for instance, = 0, — r) one can maximize the overlap over the two phases (j)'^, (j)^ . 

To solve this maximization problem it is useful to think in geometrical terms: each wave form h{f) is seen as 
a "vector" h in an infinite-dimensional vector space W endowed with the Euclidean metric defined by the (Wiener) 
scalar product. For any given t^, t^, one sees by expanding the cosines by the usual addition formula that the 
two-parameter family of "vectors" h^{C^ ,<f)f) span a 2-plane, i.e. a two-dimensional linear subspace of W. More 
explicitly, an unnormalized basis of this 2-plane is hi , ft-^ with 

ht ^ h^iC^ = = 0) , /.^ ^ (C^ = l,</,f = I) ; (B2) 

a generic vector in the 2-plane being h^{X"^) — Xf hf + /i^ with = C'^ cos(f>^, A^ = C'^ sin0^. Similarly, the 
two-parameter family of "exact" vectors can be written as {\^) — \f + with the same definitions as 

above with the label A being changed into X . Optimizing over the phases means finding the maximum over the A"* 
and A^ of 

where denotes the unit vector /i"*(A)/||/i'*(A)||. 

Directly attempting to maximize cos (A"*, A"^) is very cumbersome. The problem is, however, reduced to an 
easy one if one introduces orthonormalized bases in both 2-planes: say (ef ,e^) in the ^-(2-plane) and (e^,ef") in 
the X-one, with (e^, e^) — Sab = (ejf , e^); a,b= 1, 2. For instance, these orthonormahzed bases can be defined as 

et^\\ht\\-'h^, (B4) 

^ \\ht\r [\\htr Whir {ht h^?]-'^ mtr - (hthi) h^} , 

for the A-plane, and similarly for the A-plane. 



The overlap (B3) is then the scalar product between two unit vectors (one in each plane) which can be parametrized 
as cos 6 ai3 — i^oT^f) where — cosae^ + sinaef, — cos/3ef^ -f sin/3e^. Let Px denote the orthogonal projector 
onto the A-plane, and pa denote the orthogonal projection of e^, i.e. explicitly 

Po. = Pxiei) = {et ef > ef + (e^, e^) . (B5) 

The scalar product (e^, e^) is equal to {pa, e^)- It is maximized over /3 when is parallel to pa, in which case its 
value is the norm of pa ■ This shows that the maximum of cos 9 ax is equal to the maximum over a of the norm of Pa ■ 
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(cos Oax) 



max \\po 



(B6) 



On the other hand, 



where 



cos a pi + sm a p2 ■ 



pi ^ Px(ef ) = (ef , ef ) ef + (ef , ) , 
P2 s Px(e^) = (e^, ef ) ef + (e^, ef ) ef . 



(B7) 



(B8) 



In geometrical terms, Pa describes, as a varies, and ellipse in the X-plane (the projection of the circle e^ = 
cos a ei + sin a e2) and the maximum projection onto the X-plane corresponds to the semi-major axis. The square 
IIPqII^ = (Pq,Pq) reads 



\\Pa 



^ = A cos^ a 



B sin a + 2C cos a sin a , 



(B9) 



where 



^^|lPl|p = (ef,ef)^ + (ef,ef)^ 
B^\\p,r = {ete^? + {ete^)\ 
C ^ {p,,p,) = (ef , ef ) (ef , ef ) + (ef , ef ) (ef , ef ) . 



(BIO) 



Maximizing over a is now easy (using cos^ a = (1 + cos2q!)/2, sin^ a = (1 — cos2q!)/2, 2 sin a cos a = sin2Q? and 



maximizing over 2a) and yields finally 

(cos 9 ax) n 



A + B 



A-B 



(Bll) 



Inserting the definitions of the orthonormalized vectors Eq. (B4) into the definitions Eq. (BIC) of A, B and C, can 
allow one to express (cosf?Ax)max only in terms of various scalar products of the initial vectors /if, /if, /if, /if. It 
is easily checked that the final answer does not depend on the choice of basis in the A- and X-planes, and can (if 
wished) be expressed only in terms of the "2- forms" uo^ = /if A /if and uj-^ = /if A /if (and of the Euclidean structure 
of W). 



The result Eq. (Bll) gives the best possible overlap when optimizing separately over the phases of the exact 



and approximate signals. This gives the mathematical measure of the closeness of the two wave forms. However, in 
practice we do not have access to the phase of the exact signal. It might happen that the latter phase, i.e. equivalently 
the angle (3, is not optimum. Therefore, a physically more relevant measure of the closeness of the two wave forms 
(especially for the purpose of detection) is obtained by first optimizing over a (the parameter we can dial) and then 
considering that (3 has the worst possible value. In terms of the geometric reasoning given above one finds that the 
worst possible case corresponds to the semi-minor axis of the ellipse given by Eq. (B9), i.e. 



min max (cos Oax ) 

f} a 



A + B 



A-B 



(B12) 



In our simulations we considered both measures of th e closeness of the two signals. We use Eq. (Bll) when we study 
t he m athematical convergence and we u se Eq. ( B12 ) when we are interested in the detection. We shall refer to Eq. 
(Bll) as the best overlap and Eq. (B12D as the minimax overlap. 



[1] E. Poisson, Phys. Rev. D47, 1497 (1993). 



21 



[2] C. Cutler, L.S. Finn, E. Poisson and G.J. Sussmann, Phys. Rev. D47, 1511 (1993). 

[3] H. Tagoshi and T. Nakamura, Phys. Rev. D49, 4016 (1994). 

[4] M. Sasaki, Prog. Theor. Phys. 92, 17 (1994). 

[5] H. Tagoshi and M. Sasaki, Prog. Theor.Phys. 92, 745 (1994). 

[6] L. Blanche!, T. Damour, B.R. Iyer, CM. Will and A.G. Wiseman, Phys. Rev. Lett. 74, 3515 (1995). 
[7] L. Blanchet, T. Damour and B.R. Iyer, Phys. Rev. D51, 5360 (1995). 
[8] CM. Will and A.G. Wiseman, Phys. Rev. D54, 4813 (1996). 

[9] L. Blanchet, B.R. Iyer, CM. Will and A.G. Wiseman, Class. Quantum. Gr. 13, 575, (1996). 
[10] L. Blanchet, Phys. Rev. D54, 1417 (1996). 

[11] T.Tanaka, H. Tagoshi and M.Sasaki, Prog. Theor. Phys. 96, 1087 (1996). 
[12] C. Cutler et al., Phys. Rev. Lett. 70, 2984 (1993). 
[13] E. Poisson, Phys. Rev. D52, 5719 (1995). 
[14] E. Poisson, Phys. Rev. D55, 7980 (1997). 

[15] The most accuracy-demanding case is the wave form emitted by a binary system made of a (~ 1.4m0) neutron star and 
a (^ lOmg) black hole [0|l|]. 

[16] In the present work we shall only consider as parameters of the signal (besides an overall amplitude): two "kinematical" 
parameters tc, (a reference time and a reference phase); and two "dynamical" ones mi, m2 or equivalently m, rj. Our 
intention here is to compare the performance of new approximants versus standard ones. Evidently, our results can then 
be trivially generalized by constructing the P-approximants of more general wave forms if the need arises. 

[17] T. Damour and N. Deruelle, Phys. Lett. 87A, 81 (1981); and CR. Acad. Sci. Paris 293 (II) 537 (1981); T. Damour, 
CR. Acad. Sci. Paris 294 (II) 1355 (1982); and in Gravitational Radiation, ed. N. Deruelle and T. Piran, pp 59-144 
(North-Holland, Amsterdam, 1983). 

[18] We restrict our attention to circular(ized) systems. 

[19] We henceforth use units such that G = c = 1. 

[20] For simplicity, we consider only diagonal approximants (n, n). 

[21] Optimization over phase (p is quite subtle and it is treated at length in Appendix ^. 

[22] A loss in the overlap of 3% is equal to a reduction in the potential event rate by 10% (see Sec. 0and Ref. ^ for details). 
[23] The thresholds one may wish to retain to define the smallness of the biases will evidently depend on the context. 
[24] But not sufficient, though one will, in practice, assume it to be so. 

[25] L. Blanchet, Gravitational wave tails - of - tails, Class. Quantum Grav., submitted for publication; L. Blanchet, B. R. Iyer 
and B. Joguet, paper in preparation; L. Blanchet, G. Faye, B. R. Iyer, B. Joguet and B. Ponsot, work in progress; C. M. 
Will and M. E. Pati, work in progress; G. Schafer and P. Jaranowski, work in progress. 

[26] Note that using this rough approximation is consistent with our main purpose here which is to compare the usual post- 
Newtonian approximants of the phase evolution to new approximants. We leave to future work the problem of deciding 
with what accuracy one should correspondingly model the tensorial amplitudes of each Fourier component of the signal. 

[27] T. Damour and N. Deruelle, Ann. Inst. H. Poincare 44, 263 (1986); T. Damour and J.H. Taylor, Phys. Rev. D 45, 1840 
(1992). 

[28] The choice of v{^^ does not matter as any change in this quantity can be reabsorbed in tc and $c. However, our approximants 
to E'ji (v) will have a zero (physically corresponding to the last stable orbit) and it will be natural and convenient to choose 
v{^^ as the location of this zero. 

[29] C. Itzykson and J.B. Zuber, Quantum Field Theory (McGraw-HiU, 1980), p. 83. 

[30] L. Blanchet and T. Damour, Phys. Rev. D 37, 1410 (1988). 

[31] In other words, we are assuming some minimum "structural stability" when turning on 77. 

[32] L.E. Kidder, CM. Will and A.G. Wiseman, Class. Quantum Grav. 9, L125 (1992); Phys. Rev. D 47, 3281 (1993). 
[33] Note that we consistently label the successive approximants by the order in velocity. E.g. a 2PAi'-accurate object has the 
label 4. 

[34] N. Wex and G. Schafer, Class. Quantum Grav. 10, 2729 (1993). 

[35] It is to respect the consistency of the linearized gravitational equations (9^ T^" — 0) that we had to take explicitly into 
account the interaction accelerating the masses on their orbits. 



[36] E.g., if one factorized out the log factors in the Flux function of Eq. (4.1) as in Eq. (4.5) 

1 + in^^ lye,v' 



F{v; ri = 0) = — rj v 
5 



4- ■ 



j:f,v^ 



then the corresponding coefficients Fk are given hy Fk — Ai^ , k — ■■■ 5 ; Fk = + Bk In vi^o , fc = 6 ■ ■ ■ 11. The Ik are as 
given by Eq. (|!|) 
[37] With the convention g{f ) = J dte^'^'-f'g{t). 

[38] B.S. Sathyaprakash and S.V. Dhurandhar, Phys. Rev. D 44, 3819 (1991). 

[39] B.J. Owen, Phys. Rev. D53, 6749 (1996). 

[40] B.J. Owen and B.S. Sathyaprakash, paper in preparation. 



22 



[41] C. Cutler and E. Flannagan, Phys. Rev. D 49, 2658 (1994). 

[42] As A{h, g) is the infinite-dimensional generalization ol the cosine of the angle between two unit vectors, we could define 
the angular distance between two normalized wave forms {\\h\\ = 1 = \\g\\) by S = arccos A(h,g) and the corresponding 
chord by \\h - g\\ = 2sin | = yf2{l-A{Kg)). 

[43] Since our P4 -approximant of fiux has a pole in the region of interest we have instead used the Ps'-approximant in computing 
the overlaps and have not used the P5 -approximant (which also has a pole) at all. 

[44] Actually several equivalent distance functions such as 9 or 2 sin | as defined above ji^ . 

[45] CM. Bender and S.A. Orszag, Advanced mathematical methods for scientists and engineers (McGraw Hill, Singapore, 
1984). 

[46] The problem of phase maximisation was first considered by Schutz for the case of two wave forms which both have the 
same shape: B.F. Schutz, in The Detection of Gravitational Radiation, edited by D. Blair (Cambridge, 1989) pp 406-427, 
(also see Ref. [Q). Here we are addressing the problem of maximisation over phase when the signal and the template wave 
forms belong to different post-Newtonian family. 



23 



TABLE I. Location of the last stable circular orbit determined by the T- and P-approximants in the test mass case. The 
P-approximants predict the exact location at orders and beyond. At order the last stable orbit is not defined by 
P-approximants. 



n 


•''Iso /^Iso 






xf^^ = 0.1667 


2 
4 
6 


3.0000 
1.4415 
1.1705 




1.0000 
1.0000 



TABLE IL Overlap integrals of a test mass wave form whose Fourier transform is computed using stationary phase approx- 



imation with the same wave form but whose Fourier transform is computed using numerical fast Fourier transform, n stands 
for the order of the approximant with X denoting the exact wave form. 



n 




Ao 






Bo 


Co 






(T°,T°) 






X-'-ni-'-n/ 


/yO yO\ 

\ Tl ? -'-71 / 



2 
4 
6 
8 
10 
X 




1.000 
0.999 
0.999 
0.999 
0.999 
0.999 
0.999 






0.995 
0.994 
0.990 
0.986 
0.988 
0.987 
0.978 


0.967 
0.970 
0.955 
0.944 
0.945 
0.940 
0.918 


TA 
Eq. ( 


lBLI 
B12 


5 III. Faithfulness of the T- and P-approximants in 
)) together with the best possible overlaps (cf. Eq. ( 


the 
Bll 


test 
)) in 


mass case. Values quoted are the 
parenthesis. 


minimax overlaps (cf. 



n Ao Bo Co 

















4 


0.893 


0.868 


0.713 


0.833 


0.760 


0.964 




(0.894) 


(0.868) 


(0.720) 


(0.838) 


(0.795) 


(0.984) 


5 


0.545 


0.979 


0.338 


0.973 


0.496 


0.991 




(0.545) 


(0.980) 


(0.343) 


(0.974) 


(0.514) 


(0.998) 


6 


0.971 


0.996 


0.965 


0.988 


0.962 


0.994 




(0.972) 


(0.996) 


(0.968) 


(0.990) 


(0.968) 


(0.999) 


7 


0.989 


0.997 


0.858 


0.986 


0.891 


0.993 




(0.989) 


(0.998) 


(0.862) 


(0.988) 


(0.924) 


(0.999) 


8 


0.972 


0.998 


0.711 


0.989 


0.854 


0.992 




(0.972) 


(0.999) 


(0.717) 


(0.992) 


(0.872) 


(0.999) 


9 


0.995 


0.996 


0.928 


0.992 


0.933 


0.994 




(0.995) 


(0.996) 


(0.933) 


(0.994) 


(0.957) 


(0.999) 


10 


0.988 




0.849 




0.903 






(0.988) 


(-) 


(0.854) 


(-) 


(0.937) 


(-) 


11 


0.988 


1.000 


0.858 


1.000 


0.908 


1.000 




(0.989) 


(1.000) 


(0.862) 


(1.000) 


(0.940) 


(1.000) 
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TABLE IV. Cauchy convergence of the T- and P-approximants in the test mass case. Values quoted are the best possible 



overlaps. 



n 




Ao 






Bo 






Co 










/pO pO \ 


/^O rp(S \ 




(pO pO \ 






/pO pO \ 
X-" n> J n + 1 / 


4 
5 
6 
7 
8 
9 

10 


0.496 
0.528 
0.953 
0.988 
0.985 
0.997 
1.000 




0.918 
0.984 
1.000 
0.999 
0.999 
0.999 


0.374 
0.330 
0.770 
0.913 
0.833 
0.965 
0.997 




0.914 
0.986 
1.000 
0.997 
0.999 
0.994 


0.653 
0.541 
0.957 
0.982 
0.945 
0.985 
0.999 




0.992 
0.999 
1.000 
0.998 
0.999 
0.999 



TABLE V. Effectualness of the T- and P-approximants in the test mass case. Values quoted are the minimax overlaps 
together with the percentage bias in the estimation of total mass 100 (l — m^/m) in parenthesis. 



n 


Ao 




Bo 






Co 










{{Tlx'')) 




{{T°,X'')) 






4 0.993 




0.973 


0.971 


0.999 


0.899 




0.998 


( 0.036) 




( -0.079) 


( -18.408) 


( -3.361) 


( -2.658) 




( -0.649) 


5 0.824 




0.996 


0.986 


0.999 


0.806 




0.999 


( -0.382) 




( -0.029) 


( -104.846) 


( -1.428) 


( -15.040) 




( -0.314) 


6 0.990 




1.000 


0.983 


0.997 


0.988 




0.999 


( 0.039) 




( -0.004) 


( -2.875) 


( -0.040) 


( -0.605) 




( 0.000) 


TABLE VL Cauchy 


converf 


;ence of the T- 


and P-approximants in 


the comparable 


mass case. Values quoted 


are the best 


possible overlaps Eq. 


31^). 
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A 




B 






C 




\-'-n ■ -'n + 1 / 




X-" n 5 n + 1 / 


X^-n ) -^n+l/ 


/pv pv \ 

X-* n 5 n + 1 / 


X-* n 5 -^n + l / 




/pv pv \ 

X-* n 5 n + 1 / 


3 0.904 




0.248 


0.482 


0.181 


0.703 




0.374 


4 0.411 




0.863 


0.338 


0.872 


0.545 




0.985 


5 0.536 




0.976 


0.321 


0.979 


0.537 




0.995 


6 0.946 




0.982 


0.736 


0.996 


0.928 




0.999 
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TABLE VII. Location of the last stable circular orbit determined by the T- and P-approximants in the finite mass case for 
different values of the parameter kq- At order the last stable orbit is not defined by P-approximants. At orders v'^ and 
beyond the P-approximants predict the location of Iso pretty well. 







V 


= 1/4 








■n = 14/(11.4)2 




n 


'''Iso /-^Iso 










•^Iso /-^Iso 




•^Iso /-^Iso 


Ko = -1 








= 0.1600 








x^^ = 0.1636 




2 
4 
6 


2.8852 
1.5991 
1.2652 






1.2415 
1.0205 




2.9495 
1.5060 
1.2078 




1.0976 
1.0035 


Ko = 47/39 








= 0.1986 








= 0.1792 




2 
4 
6 


2.3240 
1.2881 
1.1020 






1.0000 
1.0000 




2.6934 
1.3752 
1.1408 




1.0023 
1.0000 


Ko =2 








= 0.2145 








xt = 0.1852 




2 

4 
6 


2.1517 

1.1925 

l.()")(icS 






0. 9258 

1. ()()12 




2.6060 
1.3306 
1.118() 




0. 9698 

1. ()()()5 
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TABLE VIII. Robustness of the T- and P-approximants in the comp arable mass case: Faithfulness. Values quoted are 
the minimax overlaps together with the best possible overlaps, Eq. (Bll), in parenthesis. System D corresponds to a binary 
consisting of stars of masses 20 Mq and 1.4 Mq . In this extreme mass ratio case the P-approximants at 0{v^) are NOT faithful 
(overlaps < 96.5%). 
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(0.999) 
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(0.996) 












= 2 
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0.557 


0.779 


0.894 


0.813 


0.896 


0.985 


0.586 


0.843 




(0.557) 


(0.780) 


(0.896) 


(0.817) 


(0.905) 


(0.989) 


(0.598) 


(0.848) 


5 


0.563 


0.989 


0.339 


0.995 


0.507 


0.978 


0.286 


0.992 




(0.563) 


(0.989) 


(0.344) 


(0.995) 


(0.531) 


(0.996) 


(0.294) 


(0.992) 


6 


0.964 


0.999 


0.883 


0.978 


0.955 


0.997 


0.939 


0.984 




(0.965) 


(0.999) 


(0.889) 


(0.978) 


(0.970) 


(0.998) 


(0.949) 


(0.992) 
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TABLE IX. Robust ness of the T- and P-approximants in the comparable mass case; Effectualness. Values quoted are the 
minimax overlaps Eq. (B12) together with the percentage bias in the estimation of total mass 100 (l — m'^/m) in parenthesis. 



71 


A 






B 


c 








{{p:i,x^)) 


((^n^^'')) 
















= -1 






4 


0.891 


0.970 


0.962 


0.997 


0.881 


0.940 




(0.268) 


(-0.089) 


(-12.550) 


(-7.582) 


(-2.217) 


(-2.779) 


5 


0.811 


0.989 


0.713 


0.996 


0.664 


0.932 




(-0.400) 


(-0.039) 


(-26.491) 


(-4.907) 


(-15.860) 


(-2.660) 


6 


0.995 


0.991 


0.982 


0.999 


0.955 


0.998 




(0.011) 


(0.025) 


(-0.869) 


(-0.104) 


(-0.999) 


(0.080) 










47/39 






4 


0.868 


0.965 


0.975 


0.999 


0.972 


0.997 




(0.214) 


(-0.111) 


(-6.961) 


(-3.618) 


(1.406) 


(-1.524) 


5 


0.831 


0.999 


0.981 


1.000 


0.810 


0.999 




(-0.429) 


(-0.032) 


(-97.241) 


(-1.118) 


(-23.654) 


(-0.205) 


6 


0.984 


0.998 


0.989 


0.999 


0.996 


1.000 




(0.025) 


(0.000) 


(1.004) 


(-0.157) 


(0.391) 


(0.019) 










= 2 






4 


0.863 


0.961 


0.968 


0.998 


0.996 


0.997 




(0.357) 


(-0.121) 


(-1.300) 


(-2.874) 


(1.699) 


(-0.450) 


5 


0.836 


1.000 


0.983 


1.000 


0.988 


0.996 




(-0.339) 


(-0.025) 


(-92.233) 


(-0.056) 


(-75.053) 


(0.000) 


6 


0.981 


0.999 


0.993 


0.995 


0.998 


1.000 




(0.032) 


(0.000) 
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FIG. 1. Schematic illustration of our methodology to compute improved templates. 
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FIG. 2. Exact energy functions (a) e{v) and (b) E'{v), in the test mass case and T- and P-approximants in the comparable 
mass (with r] = 1/4) case. Note that the comparable mass case T-approximant and P-approximant, are smooth deformations 
of the test mass function. 
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FIG. 3. Newton-normalized gravitational wave luminosity in the test particle limit: (a) T-approximants and (b) 
P-Approximants. 
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FIG. 4. Newton-normalized factored gravitational wave luminosity in the test particle limit: (a) T-approximants and (b) 
P-Approximants. 
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FIG. 5. Newton-normalized energy functions in the comparable mass case. We compare the convergence of the 
T-approximants and P-approximants. Observe that the P-approximants converge much faster to the fiducial exact energy 
than the standard approximants. 
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FIG. 6. Approximants to basic energy functions e{v) are compared with the fiducial exact energy function. Convergence of 
P-approximants is apparent. 
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FIG. 7. Newton-normalized gravitational wave luminosity in the comparable mass case. Curves are plotted for three values of 
the mass ratio (0, 14/129.96 and 1/4) with thicker curves corresponding to larger values of the mass ratio, (a) T-approximants 

and (b) P-Approximants. For comparison wc have also plotted the test particle flux. Note that the T-approximants as well as 
the P-approximants, are continuous deformations of the test mass limit. 
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FIG. 8. Histogram of fraction of events accessible using T- and P-approximants relative to the case when the true signal is 
known. 
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